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Abstract: We derive general formulas for photon and dilepton production rates from an 
arbitrary non-equilibrated medium from first principles in quantum field theory. At lowest 
order in the electromagnetic coupling constant, these relate the rates to the unequal-time 
in-medium photon polarization tensor and generalize the corresponding expressions for 
a system in thermodynamic equilibrium. We formulate the question of electromagnetic 
radiation in real time as an initial value problem and consistently describe the virtual elec- 
tromagnetic dressing of the initial state. In the limit of slowly evolving systems, we recover 
known expressions for the emission rates and work out the first correction to the static 
formulas in a systematic gradient expansion. Finally, we discuss the possible application of 
recently developed techniques in non-equilibrium quantum field theory to the problem of 
electromagnetic radiation. We argue, in particular, that the two-particle-irreducible (2PI) 
effective action formalism provides a powerful resummation scheme for the description of 
multiple scattering effects, such as the Landau-Pomeranchuk-Migdal suppression recently 
discussed in the context of equilibrium QCD. 
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1. Introduction 

1.1 Motivations and overview of previous works 

Because they interact weakly with matter, direct electromagnetic signals provide sensitive 
probes of the hot and dense matter produced in high energy heavy ion collisions Q . Once 
emitted, they mostly escape the collision zone without further interaction and, therefore, 
carry direct information about the state of the emitting system at the various stages of 
its time evolution. Recently, the WA98 collaboration has reported the first observation 
of direct photons in central Pb-Pb collisions at the CERN SPS Q. The measurement 
of direct photons and dileptons is one of the major goals of the PHENIX experiment at 
RHIC Q and is among the key observables to be studied by the ALICE collaboration at 
the LHC [|. 
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Electromagnetic probes are expected to be directly sensitive to the (partonic or hadro- 
nic) nature of the relevant degrees of freedom underlying the dynamics of the produced 
matter. They have been proposed as a promising signature of the possible formation of a 
deconfined state of matter, the so-called quark-gluon plasma ^. An intense theoretical 
activity has been devoted to study the emissive power of a quark-gluon plasma vs. that 
of a hot hadron gas in thermodynamic equilibrium (for recent reviews see 0, §)• Both 
are highly non-trivial problems, which require, in particular, a detailed understanding of 
various medium effects such as, for instance in the former case, dynamical screening ||9|, 



10|1 , or interference effects between multiple scattering processes, the so-called Landau- 



Pomeranchuk-Migdal effect |11, 12, (for a recent short review, see [|14[). To make 
contact with phenomenology, one has to model the space-time evolution of the fireball. 
For this purpose, one usually makes the simplifying assumption that, due to their multiple 
interactions, the produced partons thermalize (locally) shortly after the initial impact. The 
integration over space-time is then performed by folding the equilibrium rates, computed 
in either the deconfined or the hadronic phase, with a given hydrodynamic model (see e.g. 
p!5[]), which includes, in particular, a modelization of the equation of state in both phases 
as well as of the phase transition (for a recent review, see ]l6[). 

However, although ideal hydrodynamics is doing remarkably well in describing data 
such as elliptic flow or transverse spectra at not too high momenta [0], the assumption of 
early thermalization appears questionable, at least from the point of view of perturbative 
QCD (see e.g. [18|). Indeed, detailed theoretical studies of the early-time evolution of the 
partonic system reveal that the so-called pre-equilibrium stage may not be as short as 
usually believed. In particular, in very energetic collisions at RHIC and at LHC, where it 
is expected that the partons produced in the central rapidity region are mostly gluons, it 
appears that elementary QCD processes are not fast enough to lead to efficient chemical 
equilibration between gluon and (anti-)quark populations jl^. Moreover, recent studies of 
kinetic equilibration have shown that it is extremely difficult to reach thermal equilibrium 
on time scales lower than a few fm/c, due to the strong longitudinal expansion at very 
early times [20]. Although the issue of early thermalization is still a matter of debate 
(see e.g. |21]), it seems important to consider the influence of the non-equilibrium evolu- 
tion on the various experimental signatures of the produced matter and, in particular, on 
electromagnetic signals, which are sensitive to the earliest stages of the collision. 

There are various other motivations for studying out-of-equilibrium electromagnetic 
radiation in the context of heavy ion collisions. For instance, it is expected that a non- 
trivial semi-classical pion field configuration might form during the rapid cooling of the 
system through the chiral phase transition [^]. The electromagnetic radiation associated 
with the decay of this so-called disoriented chiral condensate has been proposed as a pos- 
sible signature of the phenomenon 24 1, which in turn might provide a useful tool to 
locate the critical point of the QCD phase diagram |25]. Another interesting motivation 
is related to the fact that, because of their sensitivity to the early stages of the collision, 
direct electromagnetic signals might bring experimental information concerning the initial 
conditions realized after the nuclear impact. 

Various aspects of out-of-equilibrium electromagnetic radiation have been addressed 
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in the literature, in particular in the context of an under-saturated partonic plasma, where 
local kinetic equilibrium is assumed to be rapidly achieved, but where quarks and gluons 
are not yet chemically equilibrated 27, In that case, the momentum distributions 
of quarks and gluons are isotropic and one can closely follow equilibrium calculations. In 
particular, one can generalize the so-called hard thermal loop (HTL) resummation and 
extend the usual analysis of dynamical screening to the non-equilibrium case [28|. So 
far, these calculations have been pushed up to the two-loop level [2^ and lead to signifi- 
cantly different results from the corresponding equilibrium calculations because of the low 
population of quarks and anti-quarks at early times The case of a system out of 

kinetic equilibrium has been much less studied. Attempts in this direction generally rely on 
folding general distribution functions, taken to be the solutions of appropriate Boltzmann 



equations, with conventional matrix elements 32]. This, however, does not provide 
a consistent description of phenomena like screening and one generally has to introduce 
an infra-red cut-off by hand. Although one can also formulate the analog of the HTL 
resummation in that case (see e.g. ||3^), it is more difficult to employ in practice because 
momentum integrations cannot be explicitly performed for arbitrary distribution functions. 

An important remark to be done concerning the above calculations is that, unlike 
their equilibrium counterpart, they do not rely on a first principle derivation of the out- 
of-equilibrium photon and dilepton production rates. For instance, the most widely used 
ansatz is to start with equilibrium-like expressions for the rates, written at a given space- 
time point, and to simply replace the Bose-Einstein and Fermi-Dirac distribution functions 
appearing e.g. in the real-time equilibrium formalism by other functions of momentum 
reflecting the out-of-equilibrium particle distribution in phase-space. This procedure would 
be exact for general stationary systems; it might also be justified when the duration of the 
emission process is short enough as compared to the typical relaxation time of the emitting 
system |^^; but it can clearly not be universally valid for computing non-equilibrium 
processes. In the general case, one expects the emission process to depend non-locally on 
the space-time history of the emitting system and the main issue is to properly formulate the 
production rate itself. An attempt in this direction has been made in Ref. ||35|] . However, the 
approach presented there is widely based on the use of standard perturbation theory, which 
is known to be badly suited for the description of non-equilibrium systems.^ Interesting 
steps towards a more satisfying formulation of the problem have been made recently in 
a slightly different context in Refs. |42, |4^, where the authors discuss possible finite 
lifetime effects for a system in equilibrium. Still, the real-time formulation of [42| suffers 
from serious physical problems such as spontaneous photon emission from the vaccum or 
infinite amounts of radiated energy [43|. These difficulties can be traced back to the issue 
of properly identifying and subtracting the unobservable virtual radiation at any finite time 
41. 



^For instance, the standard perturbative expansion is known to be poorly convergent out of equilibrium, 
due to the occurrence of spurious pinch singularities which lead, in particular, to so-called secular 

terms HJi (these singularities cancel out in thermodynamic equilibrium thanks to detailed balance 



relations 



-3- 



1.2 Summary of the present work 



In this paper, we address the question of out-of-equihbrium electromagnetic radiation from 
a very general point of view and present a first principle derivation of photon and dilep- 
ton emission rates in the framework of non-equilibrium quantum field theory. We treat 
electromagnetic interactions at lowest order in perturbation theory but, otherwise, make 
no assumption concerning the internal dynamics of the emitting system. Our approach 
generalizes the real-time formulation of Boyanovsky et al. to the case of arbitrarily 
far-from-equilibrium systems. In order to deal with the problem of virtual radiation men- 
tioned above, we propose to consistently take into account the electromagnetic dressing 
of the initial state. In particular, this allows one to unambiguously identify the virtual 
contribution at the initial time and to subtract it at once. 

We obtain general formulas relating the out-of-equilibrium production rates to the 
intrinsic dynamics of the emitting system which, at lowest order in the electromagnetic 
coupling constant, is characterized by the non-equilibrium in-medium photon polarization 
tensor. This generalizes the corresponding formulas in equilibrium [|l|, which we 

recover in the appropriate limit. In particular, we show that the description of the initial 
virtual cloud is essential in order to recover known formulas for the particular case of 
stationary systems. We also discuss the case of slowly evolving (quasi-stationary) systems 
and derive the first corrections to the static photon and dilepton production rates in a 
standard gradient expansion. 

In the last part of the paper, we make contact with recent progress concerning the 
description of far-from-equilibrium quantum fields, based e.g. on the so-called two-particle- 
irreducible (2PI) effective action formalism |^^, or on truncations of Schwinger-Dyson 



equations (for a recent review, see [49|). In particular, these methods allow for practi- 
cal, first principle calculations of the non-equilibrium dynamics in realistic particle physics 
applications |5^, We discuss how they might be employed in the present context and 
propose a formalism, based on introducing a 5-shape path in real time, which allows one 
to generalize existing methods to the case of dressed initial states. Finally, we outline 
how the relevant current-current correlator can be obtained within the 2PI effective ac- 
tion formalism. We point out that finite order approximations of the 2PI effective action 
automatically generates infinite resumations of ladder-type diagrams, which correspond to 
multiple scattering processes. Recall that the latter have recently been shown to play a 
crucial role in obtaining the full photon and dilepton production rates from an equilibrated 



QCD plasma at leading order in [12, h_ 



2. The non-equilibrium set-up 

We consider a system of interacting (bosonic and/or fermionic) fields, prepared in a given 
out-of-equilibrium initial state, and which emits photons and lepton pairs in the course 
of its time evolution. We assume that the relative energy loss is negligible so that the 
process of electromagnetic emission can be considered as a small perturbation. Our aim in 
this section is to relate, at lowest order in perturbation theory, the inclusive spectrum of 
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radiated photons and lepton pairs to the intrinsic non-equilibrium dynamics of the emitting 
system, without further approximations. 



2.1 Initial value problem: Schwinger's closed time-path 

We denote by Hs the Hamiltonian describing the internal dynamics of the emitting system 
which we assume has been prepared at time t = in an arbitrary state characterized by the 
(normalized) density matrix p^. In terms of the eigenstates |n,s) of Hg, where n denotes a 
set of relevant quantum numbers, the latter can be written 



Ps = ^Pnm\ri,s){m,s\, (2.1) 



with P^,^ = {n,s\ps\m,s) . Similarly, we denote by //gm the Hamiltonian describing free 
propagation of electromagnetic degrees of freedom, that is photons and (anti-)leptons, and 
by Pern the density matrix characterizing the initial state of the electromagnetic sector. As 
above, one can write: 

/Oem = 5^PaTI«,em)(6,em| (2.2) 

a,b 

in terms of the eigenstates |o, em) of -ffem- 

In absence of electromagnetic interactions, the most general density matrix for the 
whole system 'emitter +radiation' has the factorized form: 

p = Ps® Pcm- (2.3) 

At this order of approximation, the two subsystems evolve independently: The time- 
evolution of the emitter is driven by the interactions between its internal degrees of freedom, 
whereas the (free) electromagnetic sector undergo a trivial evolution from the initial state 
( |2.2| ). Introducing electromagnetic interactions, the total Hamiltonian reads: 

H = H, + H,^ + i^int , (2.4) 

where iJint represents the coupling of the photon field with the degrees of freedom of 
the emitting system as well as with leptons. In general, electromagnetic interactions might 
also introduce non-trivial correlations between the emitting system and the electromagnetic 
degrees of freedom in the initial state and the total density matrix p cannot be factorized 
as in Eq. (p.3|). As mentioned in the introduction, we are interested in describing the 
virtual electromagnetic cloud associated with the charged degrees of freedom of the emitting 
system. The corresponding density matrix can be formally obtained by "dressing" the bare 



initial state (2^): The dressed states |n, a) can be expressed in terms of the unperturbed 



states as follows (up to an overall phase) [^] (see also 

|n, a) = f7e(0, — oo) |n, s) ® |a, em) , (2-5) 

where the unitary "dressing" operator ?7e(0, — oo) is formally equivalent to an adiabatic 
switching-on of electromagnetic interactions from the infinite past: 

[/e(0,-cx)) = Texpj -i / due'"i?i('u)| = [/](-cx),0). (2.6) 
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t=0 



> 



Figure 1: The Schwinger closed time patli. In this formulation, the initial condition is specified 
at the initial time t = and is characterized by the dressed density matrix p. 

with e = 

Hi{t) = e^(^=+^-)* //int e-*(-f^=+-f^^-)* . (2.7) 

Using Eqs. ( |2.1[ )-( ^3|) and ( |2.5[) , one obtains the required dressed density matrix, in terms 
of the bare one: 

p = C/,(0,-oo)pC/t(0,-cx)). (2.8) 

Once the initial state has been specified, the subsequent non-equihbrium evolution 
is completely determined by the full Hamiltonian H. We shall typically be interested in 
unequal-time correlation functions of the form: 

(0(ti)---0(t„)) = Tr[pO(ti)---0(t„)], (2.9) 

where the time evolution of a given operator O is given by 

0{t) = e'^^Oe-'^K (2.10) 

For given initial conditions, specified at t = and characterized by the density matrix p, 
the above correlation functions may be computed by introducing a finite closed path in 
real time [53], as shown in Fig. ||. 



2.2 Virtual cloud vs. virtual time evolution: revisiting Keldysh's time-path 

In this subsection, we exploit the preceeding adiabatic construction, in order to replace the 
initial correlations due to electromagnetic dressing by a formally equivalent time evolution 
from an uncorrelated state of the form ( |2.3[ ) in the infinite past. This allows us to treat the 
virtual dressing of the initial state and the process of radiation during the physical time 
evolution on a similar footing and considerably simplifies the calculations.^ 

To make this more precise, we make use of the interaction picture with respect to the 
unperturbed Hamiltonian Hs + -ffem introduced above: 

Oi{t) = e*(^=+^'=™)* O e~*(^=+^'='")* . (2.11) 
Interaction picture operators are related to Heisenberg operators via: 

0{t) = U^{t,0)Oi{t)U{t,0), (2.12) 

^For the leading-order perturbative treatment presented below, it is actually sufficient that e/k <C e, 
where fc is a typical energy scale in the problem. 

^We mention that a similar procedure is employed in describing interacting fields in thermodynamic 
equilibrium. There, the non-Gaussian density matrix is dealt with by introducing a formal "time evolution" 
in the imaginary-time direction. 



-6- 



where 



U{t 



,t') = Texpj -i^^ duHi{u)^ = U^t',t). 



(2.13) 



For the sake of the argument, we focus here on the case of a two-point correlation function. 
The generahzation to higher correlation functions is straightforward. Using the form ( p.^ ) 
of the density matrix and noticing that 



we can write {t > t' > 0): 

mt)o{t')) 



Tr[p] = Tr[p] = 1 , 

Ti\p U{0, t) Oi{t) U{t, t') Oi{t') U{t', 0) 
Tr p U,{-oo,0)U{0,t)Oi{t)U{t,0)U^{0,-oo) 
xU,{-oo, 0) U{0, t') Oi{t') U{t', 0) [/e(0, -oo) 
Trip K(-oo, t) Oi{t) V,{t, t') Oi{t') V,{t\ -cx)) 



(2.14) 



[du]Hi{u)) =V}{t\t), 



where we defined the unitary operator 

V^{t,t') = Texp 
which satisfies the relation: 

v;(t,-oo) = ;7(t,o)c/e(o,-cx) 

Here, we introduced the shorthand notation: 

[du] =du[0(-n)e™ + 0(n)] 



(2.15) 



(2.16) 



(2.17) 



(2.18) 



Reading the first line of Eq. ( p. 151) from right to left, one clearly sees the time evolution 
from t = to the positive time t' , then to t > t' and, finally, back to t = 0. This is the 
essence of the closed time path formulation corresponding to the contour of Fig. |l|. Here, 
the initial state is specified at i = and is characterized by the dressed density matrix p. 
Similarly, we clearly see, from the third line of Eq. ( p. 151) , that this can be reformulated as 
a "time-evolution" along the contour of Fig. |2[ which extend to the infinite past, and where 
the "initial" state, now specified at t = -co, is characterized by the bare density matrix 
p. The closed path C along the real axis, represented on Fig. is sometimes referred to as 
the Keldysh contour ||5^ in the literature. It is useful to introduce time-ordered correlation 
functions along the contour, which can be written as (see Eq. ( p.l5| )): 



rc{0(ti)---0(t„)}) =Tr pTc[Oi{ti)---Oi{tn) j[du]Hi{u))] 



(2.19) 



where Tc represents time-ordering along C. We mention that one can equivalently formulate 
the above construction in terms of path integrals by means of standard manipulations (see 
also section ^ below). 
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t=0 



Figure 2: The electromagnetic dressing of the initial state at time t — is formally equivalent to 
a time evolution from a state characterized by the bare density matrix p at t = — oo to the physical 
initial time t = 0. This can be described by introducing the above (Keldysh) path in real time. 

In the following, we shall use standard notations for two-point correlation functions: 

= ec{t-t')C>f,{t,t')±Qcit'-t)C<f,it,t'), (2.20) 

where Oa and Oj3 are two arbitrary operators and a and /? represent a set of (Dirac, 
Lorentz, etc . . . ) indices and/or spatial variables. The plus/minus sign corresponds to the 
case of bosonic/fermionic operators respectively and the the two independent components 
are given by: 

C>^(,t,t') = {0^{t)Op{t')), (2.21) 
C<^(t,t') = {Opit')0^it)). (2.22) 

3. Out-of-equilibrium electromagnetic radiation 

We now apply the formalism described above to the calculation of the out-of-equilibrium 
photon and dilepton production rates as a function of time, at leading order in the elec- 
tromagnetic coupling constant. Assuming that no real electromagnetic radiation is present 
initially, the initial state is described by Eqs. (^) and ( p.8[ ) with 

Pern = |0,em)(0,em| . (3.1) 

The generalization of the following considerations to the case where real radiation is already 
present in the initial state is straightforward. Note that the density matrix ps describing 
the initial state of the emitting system is left arbitrary. We assume that the electromagnetic 
interaction Hamiltonian has the following form (in the Schrodinger picture): 

Hmt = J SxJ>'{x)A^{x), (3.2) 

where A^(x) is the electromagnetic field operator and 

J/^(f) = JM(f) + (3.3) 

is the total electromagnetic current, including the contributions from both the emitting 
system ( J^) and the leptonic fields if): 

f{x) = ei;{x)'yf''ip{x). (3.4) 
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Here, e denotes the electromagnetic coupling constant and 7'^ are the usual Dirac matrices. 
In the case where the emitting system contains charged scalar fields, there is an additional 
contribution - quadratic in the photon field - to the interaction Hamiltonian ( |3.2| ). We 
shall consider such a term at the end of this section. From now on, we specialize to the 
case of spatially homogeneous systems for simplicity. 

3.1 Inclusive photon spectrum at time t 

At a given time t, the single-particle inclusive photon spectrum can be expressed as:^ 

dn^/{t, k) 1 



d^xd^k (27r)3 ^ 

A — _L 



(ivW(t,^)), (3.5) 



where the sum runs over physical (transverse) polarizations. Here, 

N^^\t,k) = alit,k)ax{t,k) (3.6) 

denotes the Heisenberg picture of the photon number operator for momentum k and po- 
larization A. The corresponding creation and annihilation operators are related to the 
electromagnetic field operator and its time derivative through: 

ax{t,k) = iA{t,k) + kA{t,k) , (3.7) 



/2k 

with appropriate polarization vectors e^f^\k). Here, we defined the Fourier modes of the 
field as {V = (27r)3 ^(^^(O) is the normalization volume): 



1 



A^{t,k) = I d'xe-'''-'^A^{t,x) (3.8) 



and similarly for Afj_{t,k) = dtA^{t,k). Note that the operators a{t,k) and a^{t,k) are to 
be understood as annihilating and creating excitations above the averaged (macroscopic) 
electromagnetic field which can be non-zero in general. Therefore, the averaged value 
in Eq. ( |3.5| ) is to be understood as a connected correlator: {a) a)c = {{a) — {a^)){a — 
i^))) — (fi^ o) ~ All correlation functions we shall encounter in the following 

denote connected correlators as well and we omit the subscript 'c' for convenience. Using 
Eq. (P), we get 



E 

A=± 



where we introduced the connected photon propagator: 

G^u{t,t';k) ^ (Tc{A^it,k)Alit',k)}'j (3.10) 

^Strictly speaking, the following expression (see also Eq. ( |3.45| ) below) represents the spectrum of photons 
present in the system at time t. This can be identified with the spectrum of photons which have been emitted 
after a time t under the assumption that produced electromagnetic radiation escape the system without 
further interaction. The physically relevant limit corresponds to t ^ +00. 
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as well as the photon tensor (g^^, = diag(l, —1, —1, —1)) 

If^i^ik) = ^ ei^^*{k) e';^\k) = -g^^ - k^k^ + n^n^ . (3.11) 



For transverse polarizations, K • e'^''*^ = e~^^ = 0, where = {k^,k) is the four- 
momentum of the photon. Here, we chose the (space-like) transverse polarization vectors 
to be orthonormal: e*^^-* • e*^^'-* = —6^'^', and we introduced the time-like unit vector 
{•n? = 1) such that n ■ e^'^'> = as well as the space-like unit vector {k? = — 1) 



(3.12) 



_^ _ K^" -{K ■ n) 
^ ^{K • n)2 ' 

which satisfies k ■ e^'^^ = k ■ n = 0. For on-shell photons {K'^ = 0), one has 

" -Jitw- (3.13) 



We compute the relevant two-point function in Eq. (|3.9| ) using the method described 
in the previous section and treating electromagnetic interactions in perturbation theory. 
At each order, there are two types of contributions: those which only involve the leptonic 
current and those which involve the emitting current at least once, which we shall 
call "medium" contributions. We will show in section ^ that the former vanish identically 
when the dressing of the initial state is taken into account. Therefore, we only consider 
medium contributions here. The first non-trivial contribution to the connected correlator 
(3.10) occurs at 0{e'^) and involve the current twice. It is graphically represented on 
Fig. ^. There, the wavy lines represent the free photon propagator, that is, in Fourier 
space, the connected correlator of Fourier modes of photon field operators: 

G^^(t,t';^) = (Tc{A'i^{t,k)All{t',k)})^^, (3.14) 

where we have renamed the interaction picture photon field operators (see Eq. ( |2.11| )) such 
as to emphasize that the corresponding time evolution is driven by the free electromagnetic 
Hamiltonian alone: 

^^^(t, k) = (t, k) = e*^-^* A^'ik) e"*^--* . (3.15) 

The subscript 'em' on the bracket in Eq. (|3.14 ) indicates an averaging over electromag- 
netic degrees of freedom: (• • • )em = Tr[pem • • • ]• The blob of Fig. ^ represents the out-of- 
equilibrium in-medium photon polarization tensor: 

m''{t,t'-k) = (Tc{jt{t,k)J!:Ht',k)])^, (3.16) 

where the subscript 's' on the bracket indicates that the average value is to be computed 
with the density matrix ps of the system. Here, Js{t,k) is the interaction picture of the 
Fourier transform of the current operator (x) (see Eq. (|2.1lD ) : 

J^it, k) = Jf (t, k) = e'^^^ Jf'ik) e-^-^=* . (3.17) 
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t' 



v 



Figure 3: The in-medium contribution 5Q^^{t,t'\k) to the non-equiUbrium photon propagator 
at lowest order in the electromagnetic coupling constant e. Wavy lines represent the free vacuum 
photon propagator and the dashed blob represents the unequal-time current-current correlator 
which, at this order, corresponds to the non-equilibrium in-medium photon polarization tensor 

As in the case of the photon field above, we have renamed the interaction picture for 
this operator, which only depends on the degrees of freedom of the emitting system, to 
emphasize that its time dependence is completely determined by the Hamiltonian of the 
system. We stress, in particular, that, as far as the internal dynamics of the emitting system 
is concerned, I\.pu{t,t' ]k) represents the exact unequal-time current-current correlator. 

The lowest order in-medium contribution to the photon propagator ( p. 18 ) represented 
on Fig. I reads:^ 



6Gpu{t,t'-k) 



dudv Gfj_p{t, u; k) H'"^ {u, v; k) Gau{v, t'; k) , 



(3.18) 



with time integrals along the contour of Fig. |2|. The free vacuum photon propagator reads, 
in Feynman gauge: 

Gp,{t, t'; q) = -gp, G{t, t'; k) , (3.19) 
with (see Eqs. and (|23^ )) 



G>{t,t';k) 
G<{t,t';k) 



2k 



2k 



-ik(t-t') 



Jk(t-t') 



(3.20) 
(3.21) 



Writing explicitly the time integrals along the contour, we obtain: 



{t, t'; k) = - j [du] j [dv] Gp{t, u; k) n<^(n, v; k) Gp{v, t'; k) 

J —OO J — OO 

[du] / [dv]Gp{t,u;k)UP^{u,v;k)G<{v,t';k) 

-OO J — OO 

+ [ [du] j [dv]G<{t,u;k)Ii1^{u,v-k)Gp{v,t';k) (3.22) 



where we used the notation (2.18) and where we introduced the so-called spectral compo- 
nent of the various two-point functions: 



n%{t,t'-k) = i (n>,(i, t'- k) - n<,(i, t'- k)) 



(3.23) 



^The free part of the photon propagator does not contribute to the number of produced photons and 



simply gives a constant contribution to Eq. (3.9), which is zero in the present case (see Eq. (3.1)). Therefore, 
only the "in-medium" part to the photon propagator, Eq. (p.lq) below, contributes to Eq. (3.9). 
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and similarly for the free propagator: 

Gp{t, t'-k) = i {g> {t, t'- k) - G< {t, t'- k)) = ''"^^^"^'^ . (3.24) 

Working out the time derivatives in Eq. (^|^), we obtain the following expression for the 
single-photon inclusive spectrum at time t: 

2^1^ = /_J^^][^-]«"''^'^"''^7'^^(^)n.<.(->«;fc)- (3.25) 



We end this subsection by mentioning that the above expression (3.22) for the in- 
medium contribution to the photon propagator can be rewritten in a more compact form, 
which will prove useful for later use: 

&Q^At:t';k) = / [du][dv]\ GR{t,u;k)n<,{u,v-k)GA{v,t'-k) 

+GRit,u;k)U^^{u,v;k) G<{v,t';k) 

+G<{t, u- k) Ii%{u, V- k) Ga{v, t'; k) } , (3.26) 

where we introduced the retarded and advanced components: 

Ii%{t,t'-k) = Q{t-t')Il%{t,t'-k) (3.27) 
Ii^At,t';k) = -e(t'-t)n^,(t,t';^), (3.28) 

and similarly for the free propagator G. 

3.2 Alternative derivation and physical interpretation 

Here, we present an alternative derivation of the formula ( p. 251 ) for the inclusive photon 
spectrum, not using the virtual negative time evolution formalism introduced in the pre- 
vious section. For this purpose, we work in the basis where the density matrix of the 
system is diagonal (we use Greek letters to label the corresponding states throughout this 
subsection): 

Ps = Y,Pl\a,s){a,s\, (3.29) 

a 

where (a, s| ps |/?, s) = d^p P^. 

Let us first focus on a particular state |a,s) of this statistical ensemble and construct 
the corresponding dressed initial state \a) at time t = 0. At 0(e), the latter consists of 
a superposition of the bare zero-photon state and of all one-photon states allowed by the 
dynamics (the photon cloud). Making use of the adiabatic representation (|2.5|) , it can be 
expressed as: 



\a) 



/ dne'"i?i(ii)||a,s) (g) |0,em) 

^ J —CO ^ 

|a,s)®|0,em) + yy ^ ^ ^ $W(fc) |/3, s) ® |(A, fc), em) (3.30) 



A=0 P 
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where 

\{\,k),em) =a{{k)\0,em) (3.31) 
is the one-photon state with polarization A and momentum k. The wave-function ^^^^ (k) is 



given by the appropriate matrix element of the interaction Hamiltonian ( |3.2| ), which reads 
explicitly: 

i^fl{k) = / due^^^'^^^^{(3MJ^{uMc^.^)- (3.32) 
J-oo \ Ik 

The probability amplitude of having at least one photon with (transverse) polarization A 
and momentum k at time t is obtained by evolving the state ( p.3[)| ) up to time t with 
the full Hamiltonian H and by projecting over all possible final states containing at least 
the desired photon. At leading order in e, only one-photon final states contribute and the 
associated probability amplitude consists of two distinct pieces corresponding to the cases 
where: (1) the measured photon is actually produced during the time evolution from t = 
to t ; (2) the measured photon was already present in the initial virtual cloud. Below, we 
construct the two corresponding amplitudes at lowest order in e: 

(1) the measured photon is produced between times t = and t. This corresponds to the 
case where one starts in the zero-photon part \a, s) |0, em) of the initial state ( |3.30| ). The 
corresponding amplitude can be obtained by the following sequence: First, the unperturbed 
time evolution from the initial state |a, s) (8> |0,em) at t = to an intermediate state 
\i, s) (8) |0, em) at the intermediate time < u < t. The associated amplitude is given by the 
corresponding matrix element of the unperturbed evolution operator exp{—i[Hs + -ffem]^)) 
which reduces to: 

(i,s|e-*-^="|a,s). (3.33) 

Second, the transition from the state \i, s) (8) |0, em) to one of the possible one-photon states 
|j, s) (8> |(A,/c),em) allowed by the dynamics. The corresponding transition amplitude is 
given by the matrix element of the interaction Hamiltonian between these states: 



^ ^ > {jMJ''{k)\i,s). (3.34) 



2k 

Finally, the time-evolution from state |j, s) (8) |(A, k),em) at time u to the final state |/?, s) (8 
|(A, fc),em), at time t, with the amplitude 

where the phase factor corresponds to the free propagation of the one-photon state. 

Integrating over all possible times of emission and summing over all possible interme- 
diate states, one obtains the amplitude: 

4'L«(i'^) = 4f Tdne-^'^^*-") (/3,s|e-^^=*eW(^) • Mu,k)\a,s) . (3.36) 
v2k jo 



where we used the definition ( 3.17] ) 
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(2) the measured photon is part of the initial virtual cloud. In that case, it freely 
propagates from time t = to time t. The amphtude corresponding to the transition from 
a given one-photon states \6,s) |(A, A;), em) of the virtual cloud in ( |3.30D at time i = to 
the final state s) (8) |(A,/c),em) at time t reads: 

e-*'=*(/?,s|e-*^=*|5,s). (3.37) 

Summing over the possible states of the virtual cloud weighted with the associated wave- 



function (3.32), one obtains the final amplitude: 

^^2;Lpit,k) = ^ e-'=*(/?,s|e-^=*|5,s) x <^fjik) 
s 



1 

-j= / (ine™e-^^(*-")(/3,s|e-^^=*eW(fc)- J,(u,A;)|a,s). (3.38) 



Putting everything together we get, for the total amplitude: 

Ai\{t, k) = A[';l_,p{t, k) + 4'L/3(*, k) 
ft 



^ / [d'u]e-^'=(*-")(/3,s|e-^^=*eW(^)- Js(n,fc)|a,s). (3.39) 
2A^ J — oo 



Squaring this amplitude to get the associated probability and summing over all possible 
final states s) as well as over transverse polarizations, one obtains the desired inclusive 
spectrum corresponding to the initial state ( 3.30| ). Finally, averaging over all possible initial 



states \a,s) with the appropriate weight P^, we recover our previous result, Eq. ( p.25| ), for 
the inclusive spectrum: 

p')'^ = E^ifEEKV.^)'l (3.40) 



.X=± (3 



Therefore, we see that the various contributions to Eq. ( 3.25| ) corresponding to the 



various parts of the time integrations can be given very simple interpretations. Writing 
schematically (we omit the e-terms): 

rt j-0 j-0 i-O j-t j-t rO j-t j-t 

I [du] [dv] = du dv + du dv + du dv + du dv , (3.41) 

J ~oo J —oo J —oo J —oo Jo Jo J —oo Jo Jo 

we see, in the light of the previous discussion, that the contribution from negative times 
{u,v < 0) corresponds to the probability that the photon in the final state was already 
present in the initial virtual cloud. Similarly, the contribution from purely positive times 
(u, V >0) corresponds to the probability that the measured photon has actually been pro- 
duced between times t = and t. Finally, the cross-terms, which involve integrations over 
both negative and positive times, describe the interference between these two possibilities. 
Notice that the latter depends explicitly on time and, as a consequence, contributes to the 
actual production rate (see below). 
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3.3 Photon production rate and spectrum of produced photons 

One obtains the out-of-equiHbrium photon production rate at a given time t, by taking the 
time-derivative of Eq. (p.25|). The complete expression at lowest order in e is, therefore: 



d^xd?k {2TTf 



f [du] 2Re{e*'=(*-") ^^k) n<,(t, u;k)\ , (3.42) 



where we used the fact that: n!j^(t, i'; k) = 115/ (t', t; k). Equation ( ^.42 ) relates the photon 



production rate to the intrinsic dynamical properties of the out-of-equilibrium emitting 
system, which are characterized at this order by the unequal-time photon polarization 
tensor Ii'^^{t,u]k). This generalizes the corresponding equilibrium formula |l|, ^ (see 
Eq. dU) below). 



The inclusive spectrum of actually produced photons can be obtained from Eq. (|3.25| ) 
by subtracting the contribution from the initial virtual cloud, which is simply given by Eq. 



(|3.25|) itself, written at t = 0:^ 

dn'y{t,k) _ dn^{t,k) dn^{0,k) 
prod d^xd^k d^xd^k 



d^x d^k 



(3.43) 



This is, of course, equivalent to a direct integration of the production rate (3.42) from t = 
to t, namely: 

dn^{t,k) ^ dn^{t',k) 
d^xd^k prod Jo d'^xd^k 

Notice that, although the contribution from the virtual cloud has been subtracted, 
there still remains a time integral over infinite negative times in Eq. ( p. 42 ). As we have 



seen previously, this corresponds to interference effects with the virtual cloud. 
3.4 Inclusive spectrum of correlated lepton pairs at time t 

We now come to the case of dilepton production. The density of correlated lepton pairs 
with momenta k (lepton) and p (anti-lepton) at time t is given by 



dni+i-{t,k,p} V 



d^xd^kd^p {2Trf 



5^(iVf (t,^)<(t,pl), (3.45) 



''Strictly speaking, the dressed initial state constructed in section ^contains real radiation, which should 
also be subtracted in order to obtain the amount of produced radiation. However, the number of real 
photons in the initial state can be made negligible by choosing the parameter e such that: e/fc 2> e^, 
where A; is a typical scale in the problem. As mentioned previously, the present perturbative analysis at 
leading-order is not affected by this choice as long as e/k <C e. [As an illustration of the previous point, 
consider a time-translation invariant situation, as discussed in section ^ below. In that case, the number 
of real photon in the initial state can be shown to be ~ n^^^(w = k,k)/ke, where n^y(cj,fc) is the photon 



polarization tensor in frequency space (see Eq. (4.2)). One typically has n^^(aj = k,k) = e k f{k/k) 



where the function f{x) vanishes rapidly at large x. For instance, for a system in equilibrium at temperature 
T = k, one has f{x) ~ exp(— x). Therefore, the contribution from real radiation in the initial state can be 
neglected for fc > by keeping e/k ^ e^.] 
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where the RHS is to be understood as a connected correlator. Here V = {2tt)^ 5^^\6) is the 
quantization volume and 



N^-{t,p) = bUt,p}b^{t,p) (3.46) 
Ni\t,p) = di{t,p)d^{t,p) (3.47) 

denote the Heisenberg picture number operators for leptons and anti-leptons with mo- 
mentum p and spin projection a = zb. In terms of the Fourier modes of the lepton field 
operators 

^{t,k) = ^ I d^x(i-'^-^^{t,x) (3.48) 



V 

i){t,k) = I (fxe'^-^^{t,x), (3.49) 



'V 

the relevant correlator is given by the following connected four-point function: 

Y,{K~{t,k)N'J:{t,p)) = (tr{v;(t,fc)7o^^|^7oV'(i,fc)} 

a, a' 

f — m 



xtr{7oi^7o^(t,-p1V^(t,-;^}), (3.50) 



where m is the lepton mass, = ^/k'^ + m?, Ep = \/p'^ + m?, K>^ = {E^, k), = {Ep,p), 
If = j^K^ and = "y^P^- This correlator can be computed using the method described in 
section ||. As in the case of photons the purely leptonic contributions, involving the lepton 
current only, vanish identically. Here, we focus on medium contributions, which involve 
at least once. The first non-vanishing medium contributions occur at 0{e^) and involve 
the unequal-time current-current correlator (l3.16| ). There are two possible topologies at 
this order, one of which is represented on Fig. ^ where the lines represent the free lepton 
propagator on the contour C: 



D 



it,t';k) = (rc{V'em(t,^)V5em(t',fc)})^^. (3.51) 



As we did previously for the photon field, we have renamed the interaction picture of lepton 
field operators as follows: 

tPemit, k) = iJi{t, k) = e^^^-* ^Pik) e-*^""* , (3.52) 



and likewise for ipem{t,k). With initial conditions corresponding to Eq. ( |3.1D , the compo- 
nents of the lepton propagator are given by: 

D> {t, t'-k) = e-^^^(*-*') (3.53) 

2Ek 

D<{t,t';-k) = ^^e''''^^'-''K (3.54) 
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Figure 4: The lowest order medium contribution to dilepton production. Tlie dashed blob 
represents the O(e^) in-medium photon polarization tensor. Lepton momenta follow the fermionic 
flow. 



The contribution depicted in Fig. ^ can be expressed in terms of the lowest order 
in-medium contribution to the photon propagator dGfiu, represented on Fig. I (see Eq. 

V it,k) N'J {t,p)'j = -e^ I dudvF''{t,u,v-k,p)6g,^{v,u-q), (3.55) 

where q = k -\- p is the total momentum of the lepton pair and where 

l'"'{t,u,v-k,p) = tr|D(t,u;-pl7^D(n,t;A:)7o^^^7o 

xD{t, v; k) f D{v, t; -p) 70 7o | ■ (3.56) 

Writing explicitly the time integrations along the contour and making use of Eqs. ( ^.53|) - 
( 3.54 ), one can perform the trace over Dirac indices. A lengthy but straightforward cal- 
culation shows that all but one of the contributions on the contour vanish. A similar 
calculation shows that the other possible contribution to the fermion four-point function 
gives a vanishing contribution to dilepton production after performing the relevant projec- 
tions on positive and negative energy final states. Putting everything together, we finally 
obtain the following expression for the inclusive spectrum of correlated lepton pairs at time 
t, at leading order in the electromagnetic coupling constant: 

'mIt = Wf fj'*'^ e-'(-.«..<"-«. .-(J.p) .5„<,(„, ,1 . (3.57) 
where 

ip.{k,P) = •'^ ^ " ^i"^^ Z L_ (3.58) 

i/fc -c/p 

is the usual vacuum lepton tensor |41|. One obtains the number of correlated pairs with 
total four-momentum = {qo,Q) by integrating ( 3.57 ) over momenta with the constraints 
p -|- k = q and Ep + = qo- Using 

f cfk (fp 



(27r)3 (27r)3 



{2Tif 5(4) (Q -P-K) £^,{k,p) e^(^'^+^p)* = L^,{Q^) e^^o* , (3.59) 
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where (Q^ > 4m^) 



2^ 1 / 2m2\ / 4m^\| 



we obtain: 



d3xd4g (27r)4 



/ e-^«("--) L^^^IQ^) 5G<^{v, u; q) . (3.61) 

J — oo 



Note that the general structure of the expressions (|3]57D and (|33l|) is very similar to 



that obtained in the case of photon production, Eq. (|3.25| ). In particular, the integrand 
does not explicitly depend on t and the time-dependence only appears through the upper 
bound of the time integrations. Following the line of reasoning of section |3.2| , the various 
integrations over negative and positive times involved in Eqs. (3.57), ( p.61| ) and ( 3.22| ) 



can be given simple physical interpretations. At lowest order in e, the basic process is 
the emission of an off-shell photon which decays into a lepton pair and the probability 
amplitude for having (at least) a lepton pair in the final state can be decomposed in 
various pieces corresponding to the various possible times of emission and decay of the 
off-shell photon. For instance, the pair may already be present in the virtual cloud, which 
correspond to both times being negative, or an off-shell photon from the virtual cloud may 



decay between t = and t, etc. The various pieces of the time integrations in Eqs. ( 3.57 ), 



( |3.61 ) and ( 3.22| ) correspond either to the probabilities associated with each possibility or 



to the interferences between different possibilities. 



3.5 Dilepton production rate and spectrum of produced pairs 



As in the case of photon production, one obtains the rate of lepton pair production as a 
function of time by taking the time derivative of expressions ( 0.57 ) and ( 3.61 ) above. We 
get, for pairs of lepton and anti-lepton with momenta k and p respectively: 



Sx(Pk(Pp (27r)* 
and, for a pair with total four momentum Q: 

dn(,+e-{t,qo,q) _ 



dn]2Re{e^(^*'+^^)(*-")£^^(fc,p)Je?<^(t,u;^} , (3.62) 



f [du] 2Re|e^^o(*-") L^^(Q2) 5g<^{t^ u;q)\ . 

J —oo 



(3.63) 



The number of produced pairs at time t is obtained either by subtracting the virtual 
contribution [t = 0) from expressions ( 0.57 ) and ( 3.61| ) or by integrating the above rates 
from i = to t. For instance: 



d^x d^Q 



prod 



* dn),+i-{t',qo,q) 
d^xd^Q 



(3.64) 
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3.6 Charged scalar fields 

We end this section by considering how the above formulas for the photon and dilepton 
production rates are modified in the case where the emitting system contains charged 
scalar fields such as e.g. charged pions. This is of interest e.g. in the hadronic phase in 
heavy ion collisions, or during the out-of-equilibrium chiral phase transition, where so- 
called disoriented chiral condensates might occasionally form, leading to coherent pion 
emission |22]. The accompanied electromagnetic radiation might provide an interesting 
signature of this phenomenon |^, ^] . 

Apart from the contribution of charged scalar field to the current (^^) , the electromag- 
netic interaction Hamiltonian (|3.2|) now contains an additional contribution of the type: 



HfT = e' J d'x A^'{x)A^{x) ^\x)m , (3.65) 

where 0(x) denotes the charged scalar field operator. Repeating the above analysis for 
photon and dilepton production, it is easy to see that, at lowest order in e, the presence of 
this "two-photon" vertex simply amounts to an additional tadpole-like contribution to the 
expression ( p. 181) of the in-medium photon polarization tensor. It can be formally obtained 
from Eq. ( |3.1^ ) by replacing Il^y{u, v; q) by —i K(n) g^^y 6c{u — v), where 

/d^k 
^{4{t,k)^,{t,k)),, (3.66) 

with ips{t,k) the interaction picture of the Fourier modes of the field operator (see Eq. 

ifsit, k) = (/7i(t, k) = e*^=* ^{k) e-*^=* , (3.67) 
where {V is the normalization volume) 

ip{k) = ^ [ d^xe-'''-^(P{x). (3.68) 



Writing explicitly the time integrals along the contour, one obtains that, in presence of 
scalar fields, the tadpole contribution to Eq. ( p. 22 ) reads: 

SQU*,t';q) =g^, I [du]G p{t,u;q)K{u) G< {u,t' -q) 

tadpole 

t' 



- g^u [du]G<{t,u;q)K{u)Gp{u,t';q). (3.69) 

J —oo 

Using the expressions ( 3.20| ), ( p.21| ) and ( 3.24| ) of free propagators, it is straightforward 
to show that the tadpole term does not contribute to photon production: 

lim (dtdf+ik{dt-dt>) + k^)6g<it,t';k) =0. (3.70) 



tadpole 



However, in general it gives a non- vanishing contribution to the out-of-equilibrium dilepton 
production. We stress that this contribution results from the non-trivial time-evolution 
of the emitting system. In particular, it vanishes identically in equilibrium and, more 
generally, for any stationary state, as we will show below. It is, therefore, a genuine non- 
equilibrium effect. 



-19- 



4. Applications 



In this section, we apply the general formulas derived above to the case of time-translation 
invariant (stationary) systems. In particular, we recover known expressions for the produc- 
tion rates in thermodynamic equilibrium. We also show that a consistent description of the 
initial virtual cloud guarantees that only on-shell processes contribute to production rates 
in the case of stationary systems and that the vacuum is stable against spontaneous emis- 
sion. These provide important checks of the physical consistency of the present description. 
Finally, we consider the case of quasi-stationary systems and work out, in particular, the 
first non-trivial corrections to the local expression of the photon production rate - used 
e.g. in hydrodynamic calculations - in a standard gradient expansion. 

4.1 Recovering knovi^n formulas for stationary systems 

It is instructive to reproduce the usual equilibrium expressions for the photon and dilepton 
production rates from the general out-of-equilibrium expressions derived above. For this 
purpose, we consider the case where the emitting system undergoes a stationary evolution. 
This is the case whenever the initial condition does not break time-translation invariance, 
that is when the density matrix ps commutes with the Hamiltonian Hg (it can, therefore, 
be diagonalized in the basis of eigenstates |n,s), cf. Eq. (|2.l| )): 

In that case, all correlation functions only depend on time differences and it is useful to 
work in frequency space. In particular, one has: 

/ + 00 J 
^ l^e— (*-*)n<,(u;,fc), (4.2) 

and similarly for other two-point functions. Note the relation 

n<;(w,fc) = n<^(w,fc). (4.3) 

Plugging the definition ( |4.2| ) in Eq. ( |3.42| ) and using the identity:^ 

[du] e^°(*-") = TT 6{a) +iv(-], (4.4) 

-oo 

where V denotes the principal part, one obtains the following expression for the (time 
independent) photon production rate: 

''^k = l^-''''m.(k,k) = -^n<^(.,£). (4.5) 

In writing the second equality, we used the fact that YLfj_^ is transverse {K^^ = {k^, k)): 

K''Uf,,{ko,k)=0, (4.6) 



■^Writing F{t) = Jl^[du] e'"**""', with the definition ( |2.18| ), one finds that the time derivative dF/dt = 



hm,^o+ iee"^7(a + = 0. It foUows that F{t) = F(0) = i/{a + ie). 
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as a consequence of gauge invar iance. 

Similar considerations directly lead to the following formula for the production rate of 
dileptons with invariant mass Q (see Eq. ( ^.631) ): 



(4.7) 



d^x d^Q (27r)4 

where the function 5Q^ is defined as in Eq. ([4.2|). It can be easily expressed in terms of the 
polarization tensor by using Eq. ( ^.261) which has the form of a convolution product in 
time and which, therefore, corresponds to an algebraic product in frequency space. This is, 
of course, nothing else but the fact that time-translation invariance is equivalent to energy 
(frequency) conservation. The various components of the free photon propag ator (|3Tg| ) 
have the following expressions in frequency space: 



G<(u;,( 
G^(c^,( 
G^(c^,( 



vr ^, 

— o(u> + q 
Q 

-1 



(u + ir])^ 
-1 



(4.8) 
(4.9) 
(4.10) 



(cj — irj)^ — 

where r/ = O"*". The contributions corresponding to the last two terms of Eq. ( 3.26| ) con- 
tain G^j^{qo,q) oc 5{qo + q) and, therefore, give vanishing contributions for qo > 0. The 
additional tadpole contribution appearing in the case of scalar fields (see Eq. ( 3.69| )) van- 
ishes identically for the same reason. Using Eqs. ( 3.6C| ) and ( [4.6| ), one finally obtains 
(oem = e747r): 

2m^\ / 4m^' 



dne+i- (Q) 



Op 



1 + 



n 



{qo, 



(4.11) 



Equations. ( |4.5| ) and ( |4.11 ) are the usual expressions of the photon and dilepton produc- 
tion rates for general stationary situations. For the particular case of thermodynamic 
equilibrium at temperature T = 1/(3 these are usually expressed in terms of the spectral 
function n^,y(go,^ = 2i Imn^^(qo, ^) (see Eqs. ( p.23| ), (3.27) and ( |3.28| )) by making use of 
the detailed balance (or KMS) relation: 



^(lu{lo,q) 



eq 



e/3^«n<,(go,g) 



from which it follows that 



eq 



— Imll^ 



eq 



(^0,0 



eq 



(4.12) 



(4.13) 



Putting everything together, one recognizes the familiar expressions for the photon and 
dilepton production rates in equilibrium [|l|, 46 , 41 1 . 

As an important physical consequence, the above results guarantee that the vacuum is 
stable against spontaneous emission. Indeed, the fact that only on-shell (energy-conserving) 
processes contribute in stationary situations implies that: 



n 



{qo,q) 



0. 



(4.14) 
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which can also be obtained as the zero temperature Umit of Eq. ( [4.13| ). As a byproduct, this 
has the important imphcation that the purely leptonic contributions to the photon polar- 
ization tensor, which only involve the leptonic current and have, therefore, the structure 
of vacuum contributions, do not contribute to actual photon and dilepton production rates, 
as announced in section |^. 

4.2 Slowly evolving systems: gradient expansion and off-shell effects 

As mentioned in the introduction, a widely used ansatz for computing photon or dilepton 
production rates out of equilibrium is to use the above static formulas, Eqs ( |4.5| ) or (4.11), 



written at a given time t. One then supplement these local expressions for the production 
rates with a given modelization of the time-evolution of the system, like e.g. hydrodynamics, 
or an appropriate Boltzmann equation. This procedure might be justified if the typical 
time scale characterizing the process under study, say the formation time of the produced 
photon, is short compared to the time scale characterizing the time evolution of the system. 
In contrast, we clearly see from the expressions derived previously that in the general case 
the emission rates depend non-locally on the time history of the emitting system prior 
to the time of emission. It is therefore interesting to consider a somewhat intermediate 
situation where the system is slowly evolving (compared to the duration of the emission 
process) and to derive the first corrections to the static formulas, Eqs. (|4.5| ) and ( 4.11| ). 



With the general expressions of the previous section in hand, we are in a position to work 
out these corrections by performing a standard gradient expansion. Here, we illustrate this 
point on the case of photon production. 

Generalizing the notation introduced in Eq. ([4.2|), we introduce the Wigner transform 
in frequency space, that is a Fourier transform with respect to the time difference t — t' at 

fixed r= ^{t + t'y. 

n<.(t, t'- k)= e— n<,(T; u;, k) , (4.15) 



and similarly for other two-point functions. Note that the relations (|4.3| ) and ( [4.61) gener- 
alize: 

U<;{T;u;,k)=U<^{T;u;,k), (4.16) 

and 

K^'n^,(T■,ko,k) = 0. (4.17) 
The photon production rate ( |3.42| ) can be rewritten as: 

™ ^ £H^,(f/;..J)2Re{e-.-"..-..}. (4.1S) 

where U = ^{t + u) and where we introduced the (real) function 

'y^^^(k\ _ 

J^^iT- a;, k) = -f-Hr n< (T; k) . (4.19) 



2k (27r) 
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Writing s = t — u and U = t — ^s, and plugging the first order gradient expansion 



J^j{U ;io,k) = J^^{t; lo, k) — -s dtJ^^{t; io,k) + . . . 



(4.20) 



in Eq. (4.1g|), one obtains, after simple manipulations and making use of Eq. ([4.4|): 



dn^{t, k) 
d^xd^k 



T^{t-k,k)+V / — 



duj dtJ^^{t] u, k) 



27r {k - uj) 



(4.21) 



Using Eqs. (|]T3|), ( p?T^ ) and ( pl9| ), the first term of this expression reduces to the static 
expression ( |4.5| ) written at time t. This corresponds to the local ansatz described above. 
The second term represent the first non-trivial gradient correction arising from the finite 
duration of the emission process. As expected, it involves the off-shell part of the photon 
polarization tensor^, which in turn involves energy non-conserving elementary processes. 



5. The non-equilibrium photon polarization tensor 

In the previous sections, we have related the time-dependent photon and dilepton produc- 
tion rates to the intrinsic dynamical properties of the out-of-equilibrium emitting medium 
which, at leading order in e, are characterized by the unequal-time current-current corre- 
lator (3.16). In the last part of this paper, we focus on the latter and discuss how it might 



be computed in general non-equilibrium situations. We first consider the issue of the neg- 
ative time integration which involves the unequal-time correlation function of two current 
operators at positive and negative times. We propose a general method for the calculation 
of such unusual correlation functions, which is based on introducing a 5-shape path in 
real time. This might be combined with recently developed techniques in non-equilibrium 
quantum field theory in order to compute the explicit time dependence of the relevant 
correlation functions starting from a given initial condition. Finally, we briefly discuss how 
the photon polarization tensor can be obtained from the 2PI effective action. 

5.1 General considerations: the 5-path 

For the purpose of illustration, let us consider the case of photon production. Writing 
explicitly the various part of the time integration, the production rate ( |3.42 ) can be written 
as: 

d'wfc^ = 2Re I du e^^ft-"^ ^^(t, n; k) + j"" due™ e''^^*-") .F^(t, n; ^) | , (5.1) 

where we denoted !F^{t,u;k) = 2k{2^)'^ n^i^(t, u; k). As discussed previously, the first term 
on the RHS of Eq. (|5.l| ) describes the actual production of a photon between time t = and 
time t, while the second term corresponds to interferences with the initial virtual photon 
cloud. We see that the latter involves the correlation Il^,y{t,u; k) (see Eq. (|3.16|) ) between 



8 

expression of the tensor 7'"'(fc) has to be kept 



Notice that one cannot make use of Eq. (4.17) under the frequency integral in Eq. ( |4.2l[ ) and the full 
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two current operators at positive time t and negative time u respectively. In practice, 
starting from a given initial condition, one therefore has to compute the time evolution of 
correlation functions of the emitting system in the positive as well as in the negative time 
direction (we recall that the latter is not a physical time evolution but merely a formal way 
of treating initial state correlations due to electromagnetic dressing) . This can be done by 
introducing an appropriate contour in real time, such as, for instance, that of Fig. ^. The 
latter is, however, not well suited for practical calculation as it requires one to specify the 
initial state at t = — oo. Here we shall introduce a different contour in real time which 
allows for a more useful formulation of the problem.^ 

We consider a generic quantum field theory described by a set of fundamental (bosonic 
and/or fermionic) fields, which we collectively denote by (p. For the sake of the argument, 
we focus on the time-evolution of the system and omit all possible (Dirac, Lorentz, etc . . . ) 
indices as well as spatial variables. Let us consider a general correlation function of two 
operators O at positive and negative times. We first notice that the operator with negative 
time argument can always be brought on the right of the one with positive time argument 
in the expressions of the production rates. The relevant correlation function can therefore 
be written {t,t' > 0): 

{Os{t) 0,(-t'))s = Yl {m, s| U{0, t)OU{t, 0)UiO, -t') OUi-t', 0)|n, s) , (5.2) 

n,m 

with the time-evolution operator 

U{t,t') = e-'^^^'-''K (5.3) 

The above expression can be read as follows: Starting from a given state \n, s) at time t = 0, 
evolve the system backward in time up to time —t' and insert the operator O; Then, evolve 
the system in the positive time direction from time —t' to time t and insert O again; Finally, 
evolve the system back from time t to time t = and compute the overlap with another 
state |m, s). Repeat this for all possible states and average with the appropriate weight 
= (n, s| ps \ m,s). It is easy to show that the calculation of such correlation function 
can be formulated by means of standard field theory methods by introducing a closed 5- 
shape contour along the real time axis, such as represented on Fig. ^. In particular, this 
can be given a path-integral formulation along the 5-path, allowing one to exploit the 
machinery of functional methods, which have proven extremely powerful in recent years 
for the description of non-equilibrium systems. A general time-ordered n-point correlation 
function with either positive or negative time-arguments on the contour S can be written 
as: 

rJo.oo . . . o.fa)}) ^ l^--' ■ ■ ■ (5,4) 



^It is important to notice that the contour C introduced in section ^ only concerns insertions of elec- 
tromagnetic vertices. We are free to choose any other contour to compute the internal dynamics of the 
emitting system. 
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Figure 5: The 5-path in real time. 



where denotes time-ordering along the contour S and 0{t) = 0[(p{t)]. The classical 
action A is the time integral of the appropriate Lagrangian along the contour: 

A[f] = j dtL[ip{t)]. (5.5) 

In the above expression, J^, Dip denotes the integral over classical paths along the contour 
of Fig. ^, with fixed boundaries at t = corresponding to the field configurations cj)' and 
(j) respectively. The remaining integrals / d0d0' ((/>', s| |<A) s) represent the appropriate 
average over possible initial field configurations. Note that the above formulation along 
the 5-path automatically includes the description of "usual" correlation functions with 
only positive time arguments. Note also that because of the 5-shape of the closed time 
path, positive times are always "later" (along the path) than negative times. Therefore, 
the time-ordering Ts always bring operators with negative-time arguments on the right of 
those with positive-time arguments. 

To illustrate this formalism, let us consider the time-evolution equation for the in- 
medium propagator of the fundamental field ip:^^ 

G{x,y) = {TsWs{x)ipM})s. (5.6) 

where we have reintroduced spatial variables: x = {x^,x) and similarly for y. Here, we use 
a unified notation where ip = ip^ for bosonic fields and ip = 97^70 for fermionic fields. The 
Schwinger-Dyson equation for the propagator ( p.6| ) reads: 

G-\x,y) = Go\x,y)-J:{x,y) (5.7) 

where Ti{x,y) is the corresponding self-energy and 

is the inverse classical propagator. The latter usually takes the form of a local differential 
operator: 

iG^\x,y)=V{x)5six-y) (5.9) 



^Notice that, throughout this section, we are focusing our attention on the internal dynamics of the 



emitting system and the propagator G should not be confused with the free photon propagator (3.19) 
introduced in section ^ 



where Ssix — y) = 5s{x^ — y^) 6^^\x — y). It corresponds, for instance, to the Klein-Gordon 
operator T^lx) = —{Ox + m,'^) in the case of a scalar field with bare mass rrib, or to the 
Dirac operator T>{x) = {i(^x — M) for a fermionic field with mass M. 

As usual, one introduces the various components of two-point functions as (upper/lower 
sign corresponds to boson/fermion two-point function): 

G(x,y) = e5(xO-yO)G>(x,y)±e5(y°-xO)G<(x,y), (5.10) 

as well as the spectral component: 

Gpix,y) = i[G>{x,y) t G<{x,y)) , (5.11) 

and similarly for the self energy S(x,y). The Schwinger-Dyson equation ( ^.7| ) can be 
rewritten in the form of time-evolution equations for independent components and G"^ 
by means of standard manipulations (see e.g. [55, 51]). One obtains^^ 



V{x)G>'<{x,y) = du^p{x,u)G>'<{u,y) - du ^>'<{x,u) Gp{u,y) . (5.12) 
Jo Jo 

with the notation du = du^ J d^u. Note that these equations have the same form in 
the various cases where x^ and y^ are either both positive, both negative, or have opposite 
signs. ^■^ Moreover, we observe that the evolution equation for correlation functions with x^ 
and y^ both positive (negative) only involve unequal-time two-point functions with positive 
(negative) time arguments. In contrast, in the case where x^ and y^ have opposite signs, the 
memory integrals on the RHS of the evolution equation ( |5.12| ) involve 'positive-positive', 
'negative- negative' as well as 'positive-negative' unequal-time two-point functions. For 
a given expression of the self energy, one can, in principle, solve these time-evolution 
equations starting from given out-of-equilibrium initial conditions, e.g. along the lines of 
Refs. [55, The main advantage of the present formalism is that the initial conditions 



are to be specified at the physical initial time t = 0. One may worry about the fact that the 
negative-time integration, e.g. in Eq. (|5.l| ) above, formally extends to the infinite past. In 
practice, however, unequal-time correlation functions are damped for large time differences, 
due to non-linear scattering effects and one therefore expects the negative-time integration 
to be effectively restricted to a finite range. 

Important progress have been made in recent years in describing the far-from-equilibri- 



um dynamics of quantum fields (for a recent review, see [49[). Practicable approximation 



^'^In deriving this equation, one makes use of the equal-time (anti-)commutation relations of field op- 
erators, which imply the following relations: Gp{x,y)\^o^yO — and d^oG p{x,y)\^o^yO = S^^\x — y) for 
bosons, and Gp{x,y)\^o^yO — i'yo5'^^\x — y) for fermions. 

^^For practical calculations (see e.g. p5| it often proves useful to employ the so-called statistical two- 
point function, defined as GF{x,y) = |( {x, y) ± G^{x, y) I (and similarly for the self energy), together 



with the spectral function (5.11) as independent variables. We mention that the corresponding evolution 



equations have a similar form as Eq. (5.12), and can actually be obtained from the latter by replacing, on 
both side of the equality, G^'^ — * Gf and E^'^ Ef for the statistical function and G^'^ — » Gp and 
E^'^ Ep for the spectral function. 
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Figure 6: The complete non-equilibrium photon polarization tensor at lowest order in e in terms 
of the full in-medium two-point function (lines) and amputated four-point function F^'*-' (grey 
box) in a general theory with charged fermionic fields (for a theory with charged scalar fields there 
would be an additional tadpole- like contribution). 



schemes may be based on the so-called two-particle-irreducible (2PI) effective action for- 



maUsm [47|, or on related truncations of Schwinger-Dyson equations [^]. In particular, 
these methods have proven a powerful tool for explicit, first principle calculations of the 
time evolution of equal as well as unequal-time correlation functions. Moreover, they are 
particularly well-suited to describe the damping effects discussed above. For instance, the 
latter can be systematically taken into account in a loop [^, |5^ or ^ expan- 

sion of the 2PI effective action beyond leading-order (mean- field) approximations. These 
methods might easily be combined with the present 5-path formalism in order to study 
electromagnetic radiation in genuine non-equilibrium situations. 

5.2 The current-current correlator from the 2PI effective action 

We end this paper by discussing how the photon polarization tensor can be obtained from 
the 2PI effective action. In general, the electromagnetic current operator is a bilinear in 
the charged fields describing the emitting system and the current-current correlator ( 3.16| ) 



therefore involves the four-point function of the theory. The 2PI generating functional 
automatically provides infinite ladder-type resummations for the latter ||5^. To illustrate 
this point we consider a fermionic field theory with electromagnetic current given by (in 
the Heisenberg picture) 

Jt{x) = eMxh''Mx)- (5.13) 

The exact current-current correlator is graphically represented on Fig. |^ and can be ex- 
pressed as follows: 

U^''{x,y) = {Ts{j^{x)J^{y)})^ 

= -eHv[^>^D4x,y)YD4y,x)] - e^V'ix) : F^^) : V^y) , (5.14) 

where the first term on the RHS is the one-loop contribution, with 

D^{x,x') = (Ts[Mx)Mx')})^ (5.15) 

the full in-medium fermion propagator, and the second term is the contribution from the 



^^Here, ip and ip denotes the fundamental fields of the emitting system and should not be confused with 
the leptonic fields of section ^ 



-27- 



amputated four-point function T^^^ (see Eq. ( 5.1^ ) below), which is represented by the grey 
box m Fig. |. In Eq. (|5l4| ) above, we have defined the "vertex" function 

V^p{x;u,v) = (d^{u,x)^^'D4x,v)) (5.16) 

where a and /? denote Dirac indices, and we have introduced the notation 

V^(x) : r(4) : V^(y) ^ / dtx d^zd^; dt; VL(x; n, n) r^^^Jn, n; t;, z;) V;^.(y; t;, t;) , (5.17) 



where a sum over repeated Dirac indices is imphed.^^ Here, T^^^ denotes the connected 
four-point function amputated from its external legs: 

(Ts{Mxi)Mx2)i's{x3)Mx4)})^ = -A{xi,X2) : r(4) : A(x3,X4) , (5.18) 

where we introduced the "two-particle" function 

A{u,u;v,v) = D^{u,v) D^{v,u) (5.19) 

and where we employed a similar notation as before: 

A{xi,X2) : T^^^ : A{x^,X4) = / dndndu dv A(xi, X2; u, u) r*-^-*(n, n; v, v) A{v,v,X3, X4) . 

Js 

(5.20) 

In order to illustrate how the 2PI effective action formalism might be used to compute 
the amputated four-point function, we consider the example of a purely fermionic theory 
with classical action: 

A['ip,^p] = j d'^x(i){x)[i^ -M]V'(x) + y(V',V5)) (5.21) 

For the relevant case of a vanishing fermionic "background" field, the corresponding 2PI 
effective action F can be written as |59, 51|: 

T[D^] = -iTilnD-^ - iTiD^^D^ + F2[D^] + const , (5.22) 

where the inverse classical propagator is given by 

iDo^x, y) = {i^, - M) 5s{x - y) . (5.23) 

The exact expression for the functional F2[-D^] contains all 2PI diagrams with vertices 
described by y(V',V') aiid propagator lines associated to the full fermion propagator D^. 
The functional trace Tr includes an integration over the closed time path S as well as 
integration over spatial coordinates and summation over Dirac indices. In absence of 
external sources, the equation of motion for is obtained by minimizing the effective 
action: 



5D^{x,y) 



0, (5.24) 



^^Note that whenever Dirac indices are involved, they are ordered and summed over in the same way as 
space-time variables. In the following, we omit them for simplicity. 
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Figure 7: The Bethe-Salpeter equation for tlie four-point function T'^^^ (grey box) in a theory 
with quartic interactions. The lines represent the full propagator of the theory (D^) and the dark 
grey box corresponds to the kernel A of this integral equation (see text). 



which leads to the following Schwinger-Dyson equation for the two-point function: 

D^\x, y) = Do\x, y) - S^(a;, y; D^) , (5.25) 
with proper self-energy given by: 

S^(x,y;Z?V') = -^#7^- (5-26) 
dD^{y,x) 

Higher correlation functions can be obtained from appropriate functional derivatives 
of the 2PI effective action, which generically leads to integral equations. For instance, 
one can show that the amputated four-point function satisfies the following Bethe-Salpeter 
equation (space-time variables and Dirac indices are implicit) :^^ 

r(^) = A + A : A : F^^) (5.27) 

with kernel given by: 

A{u,v;u,v) = .p, / ■ (5-28) 



Equation ( 5.27| ) is represented graphically on Fig. ^, where the light and dark grey boxes 



represent the four-point function F^*^) and the kernel A respectively. This integral equation 
generates an infinite resummation of ladder graphs with rungs given by the kernel ( 5.2^ ) . 



Thus we see that a finite-oidev truncation of the 2PI effective action in a given expansion 
scheme automatically leads to an infinite (ladder-type) resummation for the four-point 
function through Eq. ( ^.27| ).^^ These resummations correspond to the physics of multiple 
scatterings and actually generalize the gluon ladder resummation recently employed in the 
context of equilibrium QCD at high temperature to describe the LPM effect ]l2| , ^]. To 
make this last point more transparent, let us introduce the following resummed vertex 
function (see Fig. ^-.^^ 

5s{x - u)5s{x - v) + T^^\u,v) : V^(x) , (5.29) 
in terms of which the photon polarization tensor can be rewritten (see Fig. |9|): 

n^^(x, y) = -e^ V^(x) : V^y) . (5.30) 
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+ 




Figure 8: Definition of the resummed vertex (grey circle) in terms of the bare vertex (7'') 
and the amputated four-point function T^^^ (grey square). 





Figure 9: The non-equihbrium photon polarization tensor 11^'' expressed in terms of the resummed 
vertex (grey circle). 



Using equation ( 5.27 ), it is now easy to show that the resummed vertex ( 5.29]) satisfies 
the following integral equation (space-time variables and Dirac indices are implicit) 



(5.31) 



which is graphically depicted on Fig. This, together with Eq. ( 5.30| ) for the photon 
polarization tensor is a direct generalization of the type of integral equation introduced in 
Refs. |l^ to account for the LPM effect in equilibrium QCD. One can easily convince 
oneself that the resummation of one-gluon rungs needed there would simply correspond to 
a two-loop approximation of the corresponding 2PI effective action. It is remarkable that 
the infinite set of contributions needed to obtain the full photon and dilepton production 
rates at leading order in as actually correspond to a finite, 0{as), truncation of the 2PI 
effective action. 



^^This equation plays an important role for the issue of renormalization in the 2PI formalism j6^, |6l| |. 
Similar resummations in the 2PI formalism have been recently exploited for the computation of trans- 
port coefficient in 0{N) scalar theories 

''With the present notation, we have: r''*-'(«,u) : ^(a;) = J^dudvr^^^u,v\u,v)V^{u,v-,x). 
^ Using Eq. ( 5.27 ) and the definition ( 3.16 ), we have: 



18 T 



V = 7" -f- r 



(4) 



-,(4) 



7" + A : A 



A : A : r 



(4) 



'For recent discussions of gauge theories in the context of generalized effective actions, see [^, 
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Figure 10: The integral equation satisfied by the resummed vertex (grey circle). The dark 
grey box represents the kernel A. This equation resums an infinite series of ladder graphs with 
rungs corresponding to A. 
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